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We first present a simple proof that for any bipartite lattice at half filling the RKKY interaction 
is antiferromagnetic between impurities on opposite (i.e., A and B) sublattices and is ferromagnetic 
between impurities on the same sublattices. This result is valid on all length scales. We then focus 
on the honeycomb lattice and examine the theorem in the long distance limit by performing the 
low energy calculation using Dirac electrons. To find the universal (cutoff free) result we perform 
the calculation in smooth cutoff schemes, as we show that the calculation based on a sharp cutoff 
leads to wrong results. We also find the long distance behavior of the RKKY interaction between 
"plaquette" impurities in both coherent and incoherent regimes. 



I. INTRODUCTION 

The RKKY interaction between magnetic impurities 
in a metallic environment plays a crucial role in the way 
magnetic impurities order both in the dilute and the 
Kondo lattice limits Perfect nesting of the Fermi surface 
provides a mechanism for ordering, due to the divergence 
of the spin susceptibility at the nesting wavevectori^ The 
square lattice at half filling provides the best known ex- 
ample for perfect nesting and the consequent (tt, tt) anti- 
ferromagnetic ordering of the magnetic impurities. The 
minimal form of nesting appears in the undoped single 
layer of graphite^ where the Fermi surface shrinks to two 
Dirac points. However we show that the particle-hole 
symmetry provides another mechanism for ordering of 
the magnetic impurities in a metallic environment on all 
length scales. This result can be applied to the undoped 
graphene and is of value since low energy calculations 
are only valid in the long distance limit. Furthermore, 
this result provides a test for the low energy calculations 
and we show that the calculation based on a sharp cut- 
off fails this test. We also study the RKKY interaction 
between localized spins with more complicated Kondo 
interactions. Interestingly, the behavior of the RKKY 
interaction is qualitatively different for different types of 
magnetic impurities. It will be interesting to see which 
one of the scenarios we considered - if not all - can be 
realized in the lab. 

The organization of the paper is as follows. In Sec. |TT] 
we consider the simplest Kondo perturbation, where the 
impurity is localized at a lattice site, and only has an on- 
site Kondo interaction with the conduction electron spin. 
We prove a theorem for the RKKY interaction between 
these "site impurities" in bipartite lattices with hopping 
between opposite AB sublattices at half filling. The re- 
sult is that the sign of the RKKY interaction depends 
only on whether the impurities are localized at opposite 
sublattices (antiferromagnetic) or on the same sublattices 
(ferromagnetic). The sign is dictated by particle-hole 
symmetry and is thus valid on all length scales. 

Section IIIII focuses on the half-filled honeycomb lat- 
tice with nearest neighbor hopping, and examines the 
above theorem in the long distance limit. We do a low 



energy calculation for the RKKY exchange in different 
cutoff schemes to obtain the long distance behavior of 
the exchange. As will be explained, due to the nature of 
the singularity of the spin susceptibility near the Dirac 
nodes, the use of sharp cutoff is inappropriate. We do 
the sharp cutoff calculation to demonstrate this fact. We 
then examine two smooth cutoffs to find the long distance 
behavior of the exchange. Two cutoffs are considered 
to make sure that the calculation is cutoff independent. 
They both lead to the same answer in the long distance 
limit. 

It is important to note that although both antiferro- 
magnetic and ferromagnetic exchange have the same al- 
gebraic decay power in the long distance limit, the co- 
efficient of their decays are not the same. Furthermore, 
focusing on the RKKY in a single class (i.e. the same 
or opposite sublattices) , there is also magnitude "oscilla- 
tions" of the form cos (2k/) • (R — R')) due to the non- 
analyticity of the susceptibility at the wavevector con- 
necting the two nodes. 

In Sec. IIIIBI we discuss the extension of the RKKY 
between site impurities to more general cases. In partic- 
ular, we find the qualitative behavior of the RKKY for 
magnetic impurities sitting at the center of the hexagons 
of the honeycomb lattice. We distinguish this type of 
impurity by calling it plaquette impurity (vs. site impu- 
rity) . In writing the Kondo perturbation, these plaquette 
impurities can couple coherently or incoherently with the 
conduction electrons. We find that, for incoherent Kondo 
couplings, the RKKY exchange between plaquette impu- 
rities is always antiferromagnetic. 

We then obtain the RKKY behavior for plaquette im- 
purities when they have a coherent Kondo interaction 
with the conduction electrons around the plaquette. We 
find that, due to a nontrivial phase cancelation, the 
algebraic tail, present in all other situations we consid- 
ered, vanishes. 

We conclude by summarizing our results and men- 
tioning our disagreements - specially in the sign of the 
RKKY interaction - with other attemptsii^ made in cal- 
culating the RKKY interaction in graphene. We also dis- 
cuss the direct implication of this paper for the Kondo 
lattice model on the honeycomb lattice. 



2 



II. LATTICE RESULTS 

To establish notation we quickly review the RKKY 
interactioni^ Imagine putting 2 test localized spins S 
at lattice sites i and j (i 7^ j) and assume they have a 
small on-site spin exchange interaction with the conduc- 
tion electrons spin s. This perturbs the free hopping 
Hamiltonian 

^=-^E (4cja + H.C.) (1) 
(ij>,a 

by 

5H = ]{Si-Si + SrSi) (2) 

In perturbation theory the leading interaction induced 
by this term is 

-ffRKKY = JiiSi ■ Sj, (3) 

where Jy is given by 

h^-f j dr {s7{t)s+{Q)). (4) 

To be precise Jjj is the coefficient of the . However 

due to the SU(2) flavor symmetry of the unperturbed 
Hamiltonian; J\^^ , J^-^ and Jy^ are all equal. It is 
also understood that this is the static (imaginary-time- 
averaged) part of the RKKY interaction. 

Nearest neighbor hopping on a bipartite lattice can be 
viewed as a special case of hopping between A and B 
sublattices. The particle-hole transformation 

Cia(r) ^ (-1)' Cia(r), (5) 

where c and its conjugate c are the Grassmann fields, 
leaves the Lagrangian invariant. Lattice sites are labeled 
such that i is odd on A sublattices and even on B sublat- 
tices. This transformation sends the chemical potential /i 
to —/X. However at half filling /x = and the particle-hole 
transformation is a symmetry of the partition function. 
Using Wick's theorem we have 

Jij = -f j dT (cii(r)cjx(0))(Qt(T)cjT(0)). (6) 

We drop the spin indices from now on since the Green's 
functions do not depend on the spin flavor. Since the 



ground state is particle-hole symmetric, the particle-hole 
symmetry of the Lagrangian implies 

(ci(r)cj(0)) = (-l)'+j(Q(r)cj(0)). (7) 
The above relation immediately results in 

Jii--f i-^y^'J dTG,,{T)\ (8) 

where Gi^{t) is defined to be 

Gij(r) = (ci(r)cj(0)). ^ ^ (9) 

In addition Gy (t) is real, as the space-time matrix con- 
necting the Grassmann variables is real. Therefore 

m = 

The same result is obtained if the Hamiltonian of 
Eq. ([1]) contains more general hopping terms, but only 
between AB sublattices. In summary, for any bipartite 
lattice at half filling, with hopping only between AB sub- 
lattices, the RKKY interaction is antiferromagnetic be- 
tween impurities on opposite sublattices and is ferromag- 
netic between impurities on the same sublattices. 



III. LOW ENERGY CALCULATION 

Next we focus on the honeycomb lattice with nearest 
neighbor hopping and examine the theorem we proved in 
the long distance limit. The low energy theory is known 
to be governed by a 2+1 dimensional Dirac action con- 
taining the Dirac spinors - with internal A-B flavor - 
residing near the two independent nodes ifcu. Being 
"near" the nodes is formulated by introducing a momen- 
tum cutoff and a cutoff scheme in the calculations. 

The division of the honeycomb lattice to A and B sub- 
lattices is a necessity, as the Bravais lattice has a basis 
with two sites. Let us set the distance between the near- 
est neighbor sites to be 1. In the limit \R — ^ 1, the 
RKKY interaction Jab{R — R') is calculated by Fourier 
transforming Eq. ([6]) and constraining the momentums 
to be near the Dirac nodes: 



Jab{R-R') 



D,D' 



(2^)2 



-H')gi(D- 



(11) 



Here R and R' refer to the Bravais lattice vectors. D and D' are either +1 or —1. They denote near which 
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Dirac nodes ifc^i the (space) momentums reside. A ^ 1 
is the cutoff and CA(|g|) is a function that takes care of 
cutting off the momentum. Three examples for CAd^l) 
will be given below. 

Xab il) is given diagrammatically by 



(2^ 



fc + g;D 



(12) 



We use the continuum limit "translation" of the fields 
in terms of microscopic variables provided by Eq. (45) 
and (50) of Ref. to obtain: 



xi^(0,g) 



3|q| 



64 ' 

xXb(0.9) = XaUi)^ 

1 2 

Xab (0' 9) = (fe + > 



(13) 
(14) 
(15) 
(16) 



1. sharp cutoff 
Doing the integral of Eq. pTj) for D = D' = +1 using 

A|H-H'| 



C\ results in 



^AB (-'^ ^ -^') 



1 



\R-R'\3 



dx .tVo(.t). (20) 



The Bessel function Jo (a;) is obtained by doing the angu- 
lar integration. The integrand in the above equation is 
a widely oscillating function in the limit A\R — R'\ ^ 1. 
Using the asymptotic form for Jo{x) one easily obtains 



J++(J?-i2')cx 



1 



R-R'P/ 



-^sm{A\R~R'\-7r/4). (21) 



The changed scaling form, as well as the sine oscillations, 
are generated by the sharp cutoff in momentum space. 
The sine causes sign oscillations in contradiction to the 
theorem we proved in the first section, thus making the 
sharp cutoff inappropriate. Again even without the the- 
orem; we can see that the integral is dominated by larger 
values of |q|, and that is enough to make the sharp cutoff 
inappropriate. 



In the above calculations the Dirac dispersion velocity Vc 
is set to be 1 by scaling the space dimensions. 

To calculate Jab{R — R') we examine the following 
cutoff functions : 



Cl{\q\)^e{A~\q\ 



\q\ 



-I-jI/a 



(17) 
(18) 
(19) 



where 9 is the step function. 

The issue is that the susceptibility is dominated by 
larger values of |q|, and the decay that phase integration 
causes, in the sharp cutoff scheme, is not strong enough 
to compensate that. This issue becomes more severe in 
the continuum limit A|i? — — > oo. 

Since small |q| must dominate the large distance be- 
havior, we should adopt a different cutoff scheme than 
the sharp cutoff. For the same reason, the answer should 
be universal, not just the power of decay, but also its 
coefRcient. That is why we examine two different cutoff 
schemes C\ and C\ to make sure that we find the cutoff 
independent answer. 



2. smooth cutoffs 



Because of this issue we use the smooth cutoffs given by 
Eq. m]) and HH). As expected, Cl and C^, both lead to 
the same result in the continuum limit. The few integrals 
needed for this calculation are given in the Appendix [XI 
The final result is 



Jab{R-R') 



3j2 l + cos(2ki5 • {R-R')) 



647r 



\R-R'\^ 



(22) 



where « is used since it is the continuum limit A|i2 — 
R'\ ^ oo result. We also mention the results for Xaa ■ 



The same cutoff functions C\ and C\ results in 



Jaa{R~R') 



-■f l + cos(2kD • {R~ R!)) 
64^ \R-R'\^ ' 



(23) 



(24) 



B. RKKY for plaquette impurities 



A. Examining different cutoffs 

Next we do the calculations for the three cutoff schemes 
we have considered. We first show explicitly that the 
sharp cutoff is inappropriate. We then find the univer- 
sal answer for the RKKY interaction using the smooth 
cutoffs C\ and C\. 



The perturbation we considered in Eq. ([2|) is the sim- 
plest Kondo perturbation one can consider. In general 
the localized spin could have interactions with several 
conduction electrons. In the honeycomb lattice this can 
be realized experimentally by having the localized spins 
near the center of the graphenc's hexagons. 

To study these situations, let us consider a slight 
change of notations in denoting the localized spins. The 



Greek alphabet index S" is used to denote the localized 
spins in these situations. This is in contrast to the no- 
tation Si for the on-site perturbations we started this 
paper with. To establish notation we first consider the 
simpler theoretical problem in which the impurity has an 
incoherent Kondo interaction with the conduction elec- 
tron spins around the plaquette. The perturbation to the 
free hopping Hamiltonian is then 



dH = ) T^S" -Si 



< t), 



(25) 



where the sum over a (localized spins) and i (lattice sites) 
are understood. We assume S" has interactions with 
only a few Si's. For example, if is located near the 
center of a hexagon in the honeycomb lattice, we assume 
T" has only 6 nonzero elements, corresponding to the 
sites surrounding the hexagon. 

The RKKY interaction between S"' and S^^ is 



where J"^ is given by 



ja/3 ^ j^ayPj..^ 



(26) 



(27) 



and Jij is given by Eq. ([4]). 

To be concrete let us consider the honeycomb lattice 
and imagine that the localized spins are located at the 
center of the hexagons and are widely separated. We con- 
sider the "s-wave" model in which the non-zero elements 
of T" are all 1. J°'^{R — -R') is then given by sum- 
ming over 36 terms that are given by either Eq. (|22p or 
Eq. ([M)) . Since A and B sites surrounding the hexagons 
belong to different Bravais lattice sites, the separations 
in Eq. ([22]) or Eq. ((24|) are in general given by different 
R — R' , e.g. R R ± ai, R — R' ± a2, etc. Here ai 
and a2 are the primitive vectors of the Bravais lattice . 
However in the long distance limit, all the decay factors 
can be replaced by l/|i? — i?'^. 

Treating the cosine is a bit more delicate as cos(2k£) • 
(i?— R)) is not a smooth function. It is straightforward 
to show that cos (2k£) • R) is either —1/2 or 1. If we 
decompose R into the primitive vectors 



ai = (^3/2, V3/2 
02 - (0,-Vt] , 



R = niai + na2, 



together with 



fee = I 



An 
3V3 



one finds 



cos (2k£) • R) = cos (47r(m + n) /3) , 



(28) 

(29) 
(30) 

(31) 
(32) 




FIG. 1: (Color online) The honeycomb lattice. "A" sublattice 
is denoted by odd (empty) sites. Plaquettes a, /3, and the 
sites surrounding them are labeled. Fixing i on any site and 
running j G A (around a plaquette) results in {0, 1, 2} patterns 
for Oij = mi + Wi — mj — nj (mod 3). The same set is obtained 
by scanning j G B. Thus this procedure results in the set 
{-1/2, -1/2, 1} for cos(2feD ■ -Rij) and cos(fcD ■ Rij). 



However as we explain this complication does not come 
into play and the RKKY exchange between the plaque- 
tte impurities we have considered is always antiferromag- 
netic. 

To see this one has to label the sites around the two 
hexagons in terms of (m, n) pairs. It is straightforward 
to see that the set of 



O 



(mod 3) 



(33) 



for AA and BB sublattices is the same as the set for 
opposite AB sublattices, thus antiferromagnetism pre- 
vails. It can is also found (see Fig. [T]) that these sets 
will be grouped to {0, 1, 2}. Therefore the cosine contri- 
butions of Eq. and Eq. . when summed over the 
sites around two plaquettes, vanish: to the leading order 
J"^{R- R') is given by 



r^{R-R!) 



1 



l&TT \R-R'V' 



(34) 



To summarize, the RKKY exchange for plaquette im- 
purities (in the long distance limit) is always antiferro- 



5 



magnetic. Since the center of the hexagons form a tri- 
angular lattice, the RKKY exchange for the widely sep- 
arated plaquette impurities is frustrated. 

We also mention the results for the RKKY interaction 
between a plaquette impurity and a site impurity. Again 
the cos(2k£i • {R— R')), present in the RKKY for site 
impurities, when summed over the plaquette, vanishes. 
The final result is given by 



rAR-R!) 



1 



327r \R~ R 



/|3 ■ 



(35) 



The same result is obtained when the site impurity is on 
a B sublattice. 

We finally extend the incoherent Kondo coupling of 
plaquette impurities with the conduction electron given 
by Eq. to the coherent one. The coherent Kondo per- 
turbation is more physical. It can be justified by going 
back to the origin of the Kondo model, i.e. the Ander- 
son model. In the coherent Kondo coupling the 5" has 
a Kondo interaction with a coherent sum of conduction 
electron spin. The Kondo perturbation is then given by 



6H = j 5" • 



(36) 



where Sq is a coherent sum of the conduction electron 
spins around the plaquette a, eg. s+ is given by 



(37) 



and is the set of sites surrounding the plaquette a. 
The susceptibility needed in these calculations carries 
four sublattice indices: 




(38) 



fc:D' 



since in calculating \ s^Sp y , 4 sites {i, j, k, 1} are involved. 

Here i, j G T'" and k, 1 G V^^ . We avoided using a separate 
notation for sublattice indices in xJkji (9), e.g. "i" in 

Xikji (9) should be replaced with A if i € A, etc. J"'^ is 
then given by: 



"*" is used to make the distinction between the coherent 
RKKY and the incoherent RKKY of Eq. ([M]) . 
It is straightforward to find that: 

^ ^ ^ ^iDko-{R-Rid (^D'ko-{Ri-Ri) ^ q, (40) 

where is to denote that in summing over indices, each 
site index is restricted to be in a single sublattice. This 
is just an immediate result of the observation we made 
in the previous discussions and is explained briefly in the 
caption of Fig. [1] Since 

/ (0^^^''''''"''^^^SJ^'(9O = 0,q)CM), (41) 

is independent of D and D' , Eq. (PD|) implies that the 
depcndance of J"^ vanishes: 

jfiR-R') = + O{l/\R-R'\^). (42) 

Note that the incoherent Kondo perturbation corre- 
sponds to i = j and k = 1. The sum in Eq. (PO]) . with 
this further constraint, does not factor out and it does 
not vanish [See Eq. ^] . 



IV. CONCLUSIONS 

We first proved that the particle-hole symmetry for bi- 
partite lattices determines the sign of RKKY interaction 
between site impurities on all length scales. Secondly, 
the nature of the singularity of the spin susceptibility in 
graphene invalidates the use of sharp cutoff. In the sharp 
cutoff scheme, the main contribution for the RKKY inter- 
action comes from the large momenta, thus invalidating 
the low energy theory. 

We then studied the RKKY between plaquette impu- 
rities and also between a plaquette impurity and a site 
impurity. We first considered the simpler case of having 
an incoherent Kondo perturbation. For these incoherent 
perturbations, the RKKY involving plaquette impurities 
(in the long distance limit) always ends up to be an- 
tiferromagnetic. This is coming from the fact that, in 
mediating RKKY the contributions of electrons on op- 
posite sublattices dominate the contributions from the 
same sublattices [Eq. (gH) vs. Eq. ^^]. If they had 
equal strength, the l/|i2— i2'|'^ depcndance of the RKKY 
involving plaquette impurities would have vanished. 

We then focused on plaquette impurities in the coher- 
ent Kondo interaction regime. We found that the 
algebraic tail of the RKKY, present in all other situa- 
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tions we considered, vanishes. More work needs to be 
done for finding the leading contribution for the RKKY 
in this case. 

In terms of the magnitude of the RKKY in the long dis- 
tance limit, our results can be summarized symbolically 
as 



J"^|>|J^|>|JAB|>kAA|»|J° 



a/3 I 



(43) 



There had been other attempts in calculating RKKY 
in graphene4i^ They also observed a dependance 
for the RKKY. However it was claimed that the RKKY 
is always ferromagnetic due to graphene's "semimetalic 
properties". Based on our result, this may only happen 
for the plaquette impurities with coherent Kondo inter- 
actions [See Eq. ((42|) ]. and in that case the RKKY decay 
is faster than 

Finally we mention the implications of this paper for 
the Kondo lattice model on the honeycomb lattice. We 
started the study of the Kondo-Heisenberg model on the 
honeycomb lattico^ with the hope of finding the first ex- 
amples of algebraic spin liquid (ASL) phase in the pres- 
ence of a semimetal; or a Kondo insulator-Ncel decon- 
fined quantum critical points Based on the result of 
Sec. ini Jh > will be generated in the small Jk limit of 
the Kondo lattice model 

H=-tY, {^f^'t + H-c.) +JkY.S;- Si, (44) 

and in this regard the Kondo lattice model transforms to 
the model we studied in Ref. 0. 
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APPENDIX A: INETGRAL TABLE 



(27r)3 /fc2(fc-hg)2 64|g| 





/■27r 

Jo 




(Al) 
(A2) 


lim 

a — ^oo ^ 


/ dx x^e"^/" Jo(a;) = - 


-1, 


(A3) 


lim 

a — >oo ^ 


/ dx a;^e^^/"J2(a;) = 4 


-3, 


(A4) 


lim 

a — ^oo _ 


/•oo 

/ dx x^e-''^/°'\h{x) = 
'o 


-1, 


(A5) 


lim 


/•oo 

/ dx x^e-^'^"'j2{x) = 


+3. 


(A6) 
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The ellipsis in the first equation is to denote the non- 
universal pieces of the integral. They cause either expo- 
nential decay, or faster algebraic decays, in the continuum 
limit. 
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